Ab initio prediction of magnetically dead layers in freestanding 7-Ce(lll) 
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It is well known that the surface of nonmagnetic a-Ce is magnetically ordered, i.e., 7-like. One 
then might conjecture, in agreement with previous theoretical predictions, that the 7-Ce may also 
exhibit at its surfaces even more strongly enhanced 7-like magnetic ordering. Nonetheless, our 
result shows that the (lll)-surfaces of magnetic 7-Ce are neither spin nor orbitally polarized, i.e., 
a-like. Therefore, we predict, in contrast to the nonmagnetic a-phase which tends to produce 
magnetically ordered 7-like thin layers at its free surfaces, the magnetic 7-phase has a tendency to 
form a-like dead layers. This study, which explains the suppressed (promoted) surface magnetic 
moments of 7-Ce (a-Ce), shows that how nanoscale can reverse physical properties by going from 
bulk to the surface in isostructural a- and 7-phases of cerium. We predict using our freestanding 
surface results that a typical unreactive and non-diffusive substrate can dramatically influence the 
magnetic surface of cerium thin films in contrast to most of the uncorrelated thin films and strongly 
correlated transition metals. Our result implies that magnetic surface moments of Q-Ce(lll) can 
be suddenly disappeared by increasing lattice mismatch at the interface of a typical unreactive and 
non-diffusive substrate with cerium overlayers. 



I. INTRODUCTION 



Cerium even in its bulk form has been a long-standing 
challengeable system for theorists [l], and experimen- 
talists [H, H[ due to its demonstrated variable degrees of 
4f states correlations, 0-0] resulting in various astonish- 
ing properties. [1, Q Ce has been frequently appraised 
to present or criticize strongly correlated theories, fiol— 
l2lT | The reliability of the theoretical approaches has not 
yet been fully authenticated for this simple lanthanide 
metal. [22| Cerium displays an intriguing complex phase 
diagram containing a variety of solid states. [23| In this 
diagram there are two much more prominent phases, 
i.e., isostructural a and 7 solid states. A normal Curie- 
Weiss law describes the magnetic susceptibility of the 
7-phase, while the a-phase shows an enhanced Pauli-likc 
behavior. [8l[24| The 7-Ce is a ferromagnetic system with 
localized 4f states, whereas the a-Ce is a paramagnetic 
system and its 4f states tend to be more hybridized with 
the valence bands. po| The 7-Ce crystallizes in the fee 
structure with ~ 15% larger volume than its isostructural 
a-Ce phase. The isostructural a < — > 7 phase transition 
has been intensively studied, fill I22I [2fjl [26j However, the 
physics underlying this unique transition is still under 
debate, pfjj The difference between a- and 7-phases is at- 
tributed to dissimilarity in magnitude of the hybridiza- 
tion between 4f 1 and 5d 1 6s 2 states. 0] Now let us switch 
to its surface states, which can be even more problem- 
atic than the bulk states. The surface electronic states of 
cerium can be different from that of its bulk due to the 
appearance of surface states [27| and sensitivity of cerium 
to the symmetry of its surrounding environment. 14 
Eriksson et al.[14j. within their linearized muffin-tin or- 
bitals (LMTO) calculations employing the Vosko-Wilk- 
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Nusair local density approximation (VWN-LDA), found 
a-Ce(lOO) surfaces to be 7-like, i.e., magnetically or- 
dered. They[3| from their later result anticipated that 
the surface magnetism at the 7-Ce (100) layers should be 
more enhanced 7-like than the a-Ce(lOO). The later an- 
ticipation might be a natural consequence of the Stoner's 
picture that, "the tendency toward magnetism should be 
increased near metal surfaces, because of the narrowing 
of the density of states that yields a Stoncr enhancement 
in the susceptibility" [28|. Rothman and coauthors (32) 
experimentally by growing Ce(lll) on a W(110) buffer 
confirmed the photo-emission spectra measurements of 
the Gu et al.[33j in which a-like spectra were observed 
for a lattice parameter close to that of the 7-Ce. The 
later observations [HI, [33| yield a clue to think about the 
possibility of vanishing magnetic moments at the sur- 
face of magnetic materials, the so-called dead layers p8|. 
The validity of the fact that the surface of a-Ce(lll) 
is 7-like has extensively been experimentally [27L l3fjl | and 
theoretically d, Hl| verified. Nonetheless, there are un- 
fortunately few reports on the 7-Ce thin films to justify 
about the magnetic tendency of the 7-Ce(lll) surface to 
show whether it would be 7-like or a-like. These moti- 
vated us to systematically calculate the surface magnetic 
moments of a- and 7-Ce(lll) phases including spin and 
orbital polarizations as well as spin-orbit coupling which 
can be compulsory for cerium based systems, [lij We, in 
agreement with Ref . Hil and experiment [13, show that 
the surface of nonmagnetic a-phase tends to behave as a 
magnetic 7-like Ce compound. However, our result unex- 
pectedly implies that the surface of 7-Ce(lll) is a-like. 
Indeed, our calculations demonstrate that the surface of 
magnetic 7-Ce is neither spin- nor orbital-polarized. Our 
result implies that an unreactive and non-diffusive sub- 
strate can cause a magnetic transition from magnetically 
ordered layers to a nonmagnetic layers at the surface of 
a-Ce(lll) thin films. We shall discuss how cerium can 
cause such an unexpected magnetic transition. 
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II. COMPUTATIONAL DETAILS 



This work has been carried out using the program 
package WIEN2k,[35| which allows to perform accurate 
all-electron full-potential au gme nted plane waves plus 
local orbital (APW+lo) [H, H3 band structure calcu- 
lations of solids within the density functional theory 
(DFT).[38l. l39l] The Perdew-Burke-Ernzerhof generalized 
gradient approximation (PBE-GGA)[I(| has been used 
for the exchange-correlation functional. The Ce Muffin- 
tin radii were set to Rmt = 1-8 a.u. for the surface 
calculations. The expansion of the wave functions and 
charge densities were cut off by the RMT^max = 7.0 and 
G m ax = 12 parameters, respectively. A set of 21 x 21 x 1 
(19 x 19 x 19) special k-points has been used for inte- 
grations over the Brillouin zone of the supercell (unit 
cell) in the surface (bulk) calculations. The full relax- 
ations were performed with the criterion of 1 mRy/bohr 
on the exerted forces. The relativistic effects were taken 
into account by including the spin-orbit (SO) coupling 
in a second variational procedure. Orbital polarizations 
were included to consider appropriate correlation in 4f 
Ce states using the LDA+U method. [Ill Since LDA+U 
calculations can result in different solutions depending on 
initial conditions, care has been made to use an appro- 
priate density matrix to ensure that the result has not 
been trapped in a local minimum. 



III. SUPERLATTICE SURFACE STRUCTURE 




FIG. 1: "(color online)" Visualization of slab construction: 
(a) Fee unit cell, where atoms on different planes normal to 
the [111] direction are distinguished, (b) Diagonal view of the 
fee unit cell along the [111] direction indicating (111) surfaces, 
(c) Supercell containing 5 Ce(lll) layers immersed symmet- 
rically in a vacuum space as our simulated slab. 



Isostructural a- and 7-Ce crystallize in the Fm3m 
space group symmetry with the face-centered cubic (fee) 
lattice structure, as shown in Fig. Q](a). We illustrate in 
Fig. Q] (b) a diagonal view of the fee structure along the 
[111] direction normal to the (111) planes. The z axis is 
chosen to be perpendicular on the surfaces so that the x 
and y axes are parallel to them, as shown in Fig. Q] (b) . 
Thin-films can be constructed by stacking cerium lay- 
ers along our desired direction. In order to simulate the 
free Ce(lll) surfaces, a symmetric supercell is created 
by adding vacuum spaces above and below Ce layers, as 
shown in Fig. Q] (c). To determine the vacuum thick- 
ness, total energies, work functions and exerted forces on 
the surface atoms were calculated versus various vacuum 
thicknesses. Our result shows that 25 bohr is sufficient 
for the vacuum thickness, as illustrated in Fig. [1] (c), to 
avoid interactions with the nearest neighbors of the slab 
along the Cartesian z axis. 

IV. BULK PROPERTIES 

In order to more reliably simulate the surfaces of the 
a- and 7-Ce(lll) so that the surface properties of the 
a-phase can be distinguished from that of the 7-phasc 
in consistent with experiment, we would first investi- 
gate whether the method used for the surface calcula- 



tions can properly reproduce the bulk properties of these 
two phases. The ground state equilibrium volumes were 
measured to be 28.17 A 3 for the a-Ce and 34.36 A 3 for 
the 7-Ce.[42[ The most apparent difference between these 
two phases is that the lattice parameter of the 7-phase is 
larger than that of the ct-phase. To authenticate whether 
the degree of 4f hybridization with other valence states 
can cause this change in lattice parameter, we have ap- 
plied the LDA+U method(4l| to the 7-phase and used the 
PBE-GGA[Zo| for the a-phase. The spin-polarized PBE- 
GGA including spin-orbit coupling (GGA+SP+SO) re- 
sult shows an excellent agreement between our calculated 
lattice parameter, 4.808 A, with experiment, 4.830 A.(42l] 
for the en-phase. The later obtained excellent agreement 
using solely the GGA+SP+SO, with no further band cor- 
related corrections, can be taken as an indication to the 
fact that the 4f a-Ce electrons tend to behave as band- 
like itinerant electrons in agreement with Mott transi- 
tion (MT) scenario. [HG3 The GGA+SP+SO plus cor- 
relations (GGA+SP+SO+U) among 4f Ce electrons has 
been used with literature values of 6.1 eV[l5|, |43[ and 4.4 
cV[43] for the Hubbard U parameter in the 7-phase. We 
found that the former U value, 6.1 eV, results in a com- 
pletely wrong lattice parameter for the 7-phase. For U = 
0, the minimum of the E-V curve of the localized 7-Ce co- 
incides with the one of the cv-Ce. We could calculate the 
lattice parameter of the 7-phasc using the later U value, 
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TABLE I: Spin magnetic moments inside the Muffin-tin 
sphere (MT), spin magnetic moments in the interstitial region 
(Int.), total spin magnetic moments per atom (Spin), orbital 
magnetic moment inside the Muffin-tin sphere (Orb.) and to- 
tal magnetic moments per atom (Tot.) in \ib versus number 
of layers (N) within the GGA+SP+SO for the a-Ce(lll) thin 
films. We calculated total spin moment as Spin = MT+Int., 
for N = I and as Spin = MT+|lnt., for N>2. Total magnetic 
moment is calculated to be Tot. = Spin + Orb., for every N. 
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FIG. 2: "(color online)" Work function, <j) ( e V), versus 
number of cerium layers for a-Ce(lll), left vertical axis, 
and 7-Ce(lll), right vertical axis, within GGA+SP and 
GGA+SP+SO+U. 



which was corresponded [13 to the selected atomic config- 
urations of U = E{fS°)+E{f 1 S' 1 )-2E{f' 1 S 1 ) = 4.4 eV 
for the promoted localized 4f electrons to more delocal- 
ized 6s states. Our calculated lattice parameter of 7-Cc, 
5.169 A, using U = 4.4 eV is in excellent agreement with 
experiment, 5.160 A, as well. [42j We obtained 11 kbar for 
the a O 7-phasc transition pressure which is comparable 
with experimental value of 8 kbar and other theoretical 
results, [a. [l3l l44l |45| Total magnetic moments were cal- 
culated to be 0.000 and 1.142 \ib for the a- and 7-phases, 
respectively, in agreement with the SIC-LSD results. [44l| 



V. SURFACE PROPERTIES 
A. Work Function 

Insofar as the accuracy of our surface calculations are 
concerned, our success in reproducing bulk properties 
may not seem to be sufficient. Thus, we would second 
offer an estimation concerning the accuracy of our sur- 
face calculations by reporting on the work function as 
an extremely surface sensitive quantity [46j|. Therefore, 
we have also calculated the work function ,cj>(eV), as the 
minimum energy required to liberate an electron from the 
Fermi level (Ep) to a point with negligible kinetic energy 
at the center of the vacuum of the slab. The work func- 
tions were obtained using the relation of <f> = E vac — Ep, 
where E vac is estimated by the averaged electrostatic 
Coulomb potential at the midpoint of the vacuum of the 
slab and Ep is the corresponding Fermi energy. 

The work functions are given in Fig. [5] as a function 



of the film thickness for both of the a- and 7-C(lll) 
clean surfaces. As shown in Fig. [21 the SO+U causes to 
reduce the work function for both of the phases. The 
result shows that the work functions of the 7-phase is 
well converged through our five thin layers. The conver- 
gence of the 7-phase, as shown in Fig. [2j is more reliable 
than that of the a-phase. Fluctuation of the a-phase be- 
tween fourth and fifth layers shows also a small change 
of about 2% in the work functions. The result, Fig. [21 
indicates that the fluctuation of the work function is not 
substantially affected by the spin-orbit coupling and/or 
LDA+U correlations. The result shows that the work 
function oscillate as a function of the number of layers 
with the period of A = 2 layers. The oscillation indicates 
how the quantum size effects (QSEs) [47| can affect the 
results. The period of oscillation, as shown in Fig. [21 
is less than the period of the "abcabc" pattern of the 
fee stacking. This can be taken as an indication to the 
fact that here oscillations of the work function originates 
more dominantly from the quantum size effects (QSEs) 
at nanoscale than the period of thin films stacking due to 
the symmetry of the structure |47|. Our result shows that 
the work functions of the 7-phase are larger than those 
of the a-phase, for every number of layers. The later re- 
sult is in agreement with the LMTO result, [lH where the 
work functions were calculated to be 3.5 eV for a-Cc(100) 
and 4.2 eV for 7-Ce(100). Our calculated work functions 
are closer to the experimental value, 2.9±0.2,[48| than 
the LMTO values. These results, in accord with the 
MT[H G3 and Kondo volume collapse (KVC)[T3 pic- 
tures, which reconfirm that the 4f electrons tend to be 
more localized in the 7-Ce than that of a-Ce, may ensure 
that our subsequent surface results might be reliable as 
well. 
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B. Magnetic Moment 

Turning to the goal of this paper, we first report on the 
magnetic moments for the a— Ce(lll) thin films and then 
for the 7— Ce(lll) thin films. The magnetic moments are 
listed in Tab.Ufor the a— Ce(lll) thin films. The result 
shows that the direction of the total spin moments is op- 
posite to that of the orbital moments. However, total 
magnetic moments are still considerable for the surface 
of a-Ce. This implies that the surface of the nonmag- 
netic a-Ce is 7- like in agreement with experiments p7H30j 
and other theoretical results 0, Hl[ as well as the pi- 
oneer work of Eriksson and coworkers fl4). The result, 
as presented in Tab. U shows that the total magnetic 
moment changes by only a tiny percentage of about 1% 
on going from fourth to the fifth layer. Here in fact, 
cancellation errors between spin and orbital contribu- 
tions cause such a tiny percentage in the total magnetic 
moments. The surface of a-Ce(lll) is well known to 
be 7-like experimentally [27l. l30j and theoretically l3lj|. 
Therefore, the result presented in Tab. U can qualitatively 
corroborate the idea that the surface of a-Ce(lllVis mag- 
netically ordered in agreement with theoryj^, l3l| and 
experiment [27l, [30l] . For the 7-Ce(lll) thin films, we did 
not present similar tabic, because our result shows that 
such a table contains nothing more than zero moments. 
The accuracy to which we have calculated the magnetic 
moments is ±0.01 \ib for the 7-Cc(lll). Thereby if wc 
present a table for the magnetic surfaces of 7-Ce(lll), 
the table contains ignorable values of 0.00 ± 0.01 fiB- 
Such ignorable values are obtained for all components of 
the total magnetic moments, i.e. according to the abbre- 
viations given and defined in the caption of Tab. HI MT, 
Int., Spin, Orb. and as a result Tot., from N = 1 to N = 5. 
Therefore in essence, inversely the surface of the magnetic 
7-Ce is found within our ab initio DFT calculations to be 
a-like. Thus, we predict that the free surfaces of the 7- 
phase constitute dead layers at zero temperature. Such a 
contradictory situation in these isostructural phases con- 
cerning different magnetic behaviors of their bulks when 
compared with their free (111) surfaces can be consid- 
ered as another physical property for cerium. Our result 
confirms the experimental works pre formed by Rothman 
and coauthors [32j and Gu et al.[33| where a-like spec- 
tra were measured for a lattice parameter close to that 
of the 7-Ce. This represents another success of the den- 
sity functional theory (DFT) in fundamentally predicting 
complicated 4f-electron nanosystems from first principles 
of quantum mechanics without assuming any experimen- 
tal data. It is to some extent hard to expect that there 
can be found another case in nature similar to cerium 
with two isostructural phases and opposite magnetic sur- 
faces compared to their corresponded bulk counterparts. 
This implies that cerium due to its demonstrated capac- 
ity in reproducing two entirely opposite magnetic surface 
states might be considered as a unique case. Thereby, Ce 
once more serves as an interesting case to appraise the 
validity of our fundamental understanding based on the 



density functional theory (DFT). 

VI. EFFECT OF LATTICE MISMATCH 

From our hypothetical freestanding surface study one 
may also deduce that a lattice mismatch at the interface 
with a typical unreactive and non-diffusive substrate can 
play an important role on the magnitude of the magnetic 
moments in Ce thin films. The sensitivity of cerium thin 
films to the effects of lattice mismatch can be inferred by 
the fact that the differences between a- and 7-phases may 
be resulted from the larger lattice parameter of the 7-Ce 
than that of the a-Ce. The unreactive and non-diffusive 
substrate affects the lattice mismatch of the cerium over- 
layers. The lattice mismatch at the interface of a typical 
Ce(lll)/substrate causes to change lattice parameter of 
the cerium thin films. The change of lattice parameter 
in Ce can give rise to transition from a- to 7-phase. In 
uncorrelatcd thin films the change of lattice mismatch 
to some extent may not drastically change the electronic 
structures. This could be also the case for strongly cor- 
related materials. For instance, the effects of gold and 
cooper substrates on the ferromagnetic transition metals 
were refuted in leading to magnetically dead layers. (49j 
Lattice mismatch on cerium is important not only due to 
its 4f-electrons, but also more importantly due to exist- 
ing isostructural a- and 7-phascs in the vicinity of each 
other. This point might be of significant importance in 
nanotcchnology, as in this case one can control the surface 
magnetic moments ranging from enhanced magnetically 
ordered surfaces to magnetically dead layers by changing 
unreactive and non-diffusive substrates. Therefore, we 
predict that more lattice mismatch in a-Ce thin films can 
give rise to less magnetic surface. The magnetic surface 
of a-Ce thin films can suddenly disappear by increas- 
ing lattice mismatch towards lattice parameter of 7-Cc. 
The later point shows that how an unreactive and non- 
diffusive substrate can significantly influence the cerium 
thin films. 



VII. DENSITY OF STATES (DOS) 

To realize the underlying physics of the above dis- 
cussed property, we give in Figs. [3J the density of states 
(DOS) curves. As shown in Fig. [3J (a), the up and down 
bulk DOSs are symmetric with respect to each other for 
the a-Ce, which is not the case for the 7-Ce, as shown 
in Fig. [3J (f). Therefore, one anticipates to find (zero) 
nonzero magnetic moment for the (a-phase) 7-phase. As 
can be seen from Figs. [3J (a) and (f), including correla- 
tions among 4f-electrons cannot affect the later anticipa- 
tion. The calculated surface DOS curves within LDA are 
shown in Figs. [3J (b) and (g) for the a- and 7-Ce(lll), 
respectively. One notices from the later surface DOSs 
that the a- and 7-phases would exchange their magnetic 
behaviors by going from their bulks to their surfaces. As 
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FIG. 3: "(color online)" Bulk and free surface DOSs of the 
a-Ce (left panel) and 7-Ce (right panel). 



shown in Fig. [3J (g), the up and down bulk DOSs are 
symmetric with respect to each other for the 7-Ce(lll) 
surface, which is not the case for the a-Ce(lll), as shown 
in Fig. |3] (b). Therefore, in contrast to what has been 
already anticipated above for the bulk counterparts, at 
the surface the a-Ce(lll) is magnetic whereas 7-Ce(lll) 
is nonmagnetic. This shows that the above mentioned 
anticipation can be affected by the surface states. This 
means that the surface of the a-Ce is 7- like whereas 
the surface of the 7-Ce is a-like. As can be seen from 
Figs. [3] (c) and (h), including spin-orbit (SO) coupling 
and correlations among 4f-clcctrons cannot change back 



the reversed surface magnetic properties to that of the 
bulks. To shed light into the electronic structures of inte- 
rior layers, DOSs at various layers of our slab containing 
totaly 5 symmetric layers were also calculated within the 
GGA+SP+SO+U and the results are given in Figs.[3j(d) 
and (i) for the a- and 7-Cc(lll), respectively. The result 
shows that all the layers, i.e., deeper and/or shallower 
layers, arc magnetic for the a-Cc(lll) and nonmagnetic 
for the 7-Ce(Ill). The calculated Partial DOSs employ- 
ing GGA+SP+SO+U are shown in Figs. [3J (e) and (j) 
for the a- and 7-Ce(lll), respectively. The result, as 
can be seen from Fig. [3J (e), shows that the surface mag- 
netic properties of the nonmagnetic a-phase is due to the 
partial 4f-states. The result, as shown in Fig. |3J(j), also 
clarify the role of 4f-states in suppressing the magnetic 
moments at the surface of magnetic 7-phasc. 



VIII. CONCLUSION 

In conclusion, we have presented an evidence in fa- 
vor of entirely different magnetic orderings at the a— 
and 7— Ce(lll) clean surfaces as opposed to their bulk 
counterparts. Bulk magnetic moments were calculated 
to be zero and 1.14 fj.B for isostructural enhanced Pauli- 
like a— Ce and Curie- Weiss 7— Ce phases, respectively, in 
agreement with other theoretical calculations. Inversely, 
surface magnetic moments were calculated to be nonzero 
for the a— Ce(lll) and zero for the 7— Ce(lll) thin film 
layers. The contradictory situation that arises from dif- 
ferent magnetic transitions in these isostructural phases 
from their bulks to their free (111) surfaces can be con- 
sidered as another unexpected property for cerium and 
attributed to different 4f hybridizations with the other 
valence states at diverse surface circumstances. Our re- 
sult predicts that the nonmagnetic a-phasc tends to pro- 
duce magnetically ordered 7-like surface layers, whereas 
the magnetic 7-phase has a tendency to form a-like dead 
layers. In summery, we have shown that the surface of 
a-Ce is 7-like, while the surface of 7-Ce is a-like. Within 
the later result one concludes that how can cerium be of 
significant importance in nanotechnology, because phys- 
ical properties of cerium thin films can be drastically in- 
fluenced by those of unreactive and non-diffusive sub- 
strates which can impose large lattice mismatch. Our re- 
sult implies that the significant effects of unreactive and 
non-diffusive substrates on the adatoms can be expected 
solely from cerium due to its 4f-electrons and isostruc- 
tural a- and 7-phases as well as the discussed entirely 
different magnetic ordering of the isostructural phases at 
their corresponded surfaces. 
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